Quantum computing with nearest neighbor interactions and error rates over 1% 
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Large-scale quantum computation will only be achieved if experimentally implementable quantum 
error correction procedures are devised that can tolerate experimentally achievable error rates. We 
describe a quantum error correction procedure that requires only a 2-D square lattice of qubits 
that can interact with their nearest neighbors, yet can tolerate quantum gate error rates over 1%. 
The precise maximum tolerable error rate depends on the error model, and we calculate values in 
the range 1.1-1.4% for various physically reasonable models. Even the lowest value represents the 
highest threshold error rate calculated to date in a geometrically constrained setting, and a 50% 
improvement over the previous record. 



Building a quantum computer is a daunting task. En- 
gineering the ability to interact nonlocal pairs of qubits is 
particularly challenging. All existing quantum error cor- 
rection (QEC) schemes capable of tolerating error rates 
above 1% assume the ability to deterministically interact 
pairs of qubits separated by arbitrary distances with no 
time or error rate penalty [l|-|3j . The most recent of these 
works estimates a threshold error rate p t h of 5% 

It is far more physically reasonable to assume a 2-D 
lattice of qubits with only nearest neighbor interactions, 
proposed realizations of which exist for ion traps jij , op- 
tical lattices [f| , superconducting qubits @ , optically ad- 
dressed quantum dots @, S| s NV centers in diamond Q 
and many other systems. For such proposals, the leading 
QEC scheme 1(J III, which is based on the Kitaev sur- 
face code 
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12j . has been shown to possess a p t h of 0.75% 



We increase this to 1.1-1.4%, depending on 
the error model, bringing the threshold for geometrically 
constrained quantum computing above 1% for the first 
time. We achieve this by carefully using the given error 
model to calculate approximate probabilities of different 
error events and removing the need to initialize qubits. 

To overview, we will begin by describing stabilizers and 
our simplified quantum gate sequence, followed by a de- 
tailed discussion of how probable different error events 
are and how this information can be fed into the classi- 
cal decoding algorithm. We then present the results of 
detailed simulations, which apply two-qubit depolarizing 
noise with probability p2 after two-qubit quantum gates, 
single-qubit depolarizing noise with probability pi after 
identity gates, and make use of measurement gates that 
report and project into the wrong eigenstate with prob- 
ability pm . In addition to the standard error model with 
P2 = Pi = Pm = P, which we focus on, we simulate a 
balanced error model with pj = 4^2/5 and pu = %Pi/3, 
ensuring idle qubits have the same probability of error 
as a single qubit involved in a two-qubit gate and taking 
into account the fact that a measurement is only sensi- 
tive to errors in one basis. We also simulate the case 
Pi = P2/IOOO and pm = P2/IOO, modeling typical error 
ratios in an ion trap. 

A stabilizer [lH of a state \H>) is an operator S such 
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FIG. 1: a) 2-D lattice of data qubits (circles) and syndrome 
qubits (dots) and examples of the data qubit stabilizers, b) 
Sequence of CNOTs permitting simultaneous measurement of 
all stabilizers. Numbers indicate the relative timing of gates. 
The highlighted gates can be tiled to fill the plain. 
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FIG. 2: Circuit determining the sign of a stabilizer a) 
XXXX, b) ZZZZ, without explicit initialization gates. 



that S \ ^>) = \*&). An error E that anticommutes with S 
can be detected as SE |#) = -ES |\&) = —E \ty). Exam- 
ples of surface code stabilizers [l2[ are shown in fig. [TJi. 
Circuits measuring these stabilizers without explicit ini- 
tialization gates are shown in fig. [2] We assume quan- 
tum nondemolition measurements, which have been ex- 
perimentally demonstrated using ion traps optical 
lattices [l?} , superconducting qubits 18 1 and NV centers 
in diamond [19| and theoretic ally proposed for optically 
addressed quantum dots |2(| . The initial and final 
measurements match when +E is measured and differ 
when — E is measured. An appropriate sequence of two- 
qubit gates for measuring all stabilizers across the lattice 
simultaneously is shown in fig. [T]d. Data qubits execute 
identity gates while the syndrome qubits are measured. 
Repeatedly executing the gates of fig. HJ), along with 
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appropriate syndrome qubit measurements and data 
qubit identity gates, generates points in space and time 
where the stabilizer changes sign, indicating local errors. 
Renormalization techniques exist capable of processing 
perfect syndrome information 22 , 23| , however at present 
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only the minimum weight perfect matching algorithm 
2J, |25[ can be used to process the output of realistic 



quantum circuits. 

The minimum weight perfect matching algorithm takes 
coordinates and a measure of separation and matches 
pairs of coordinates such that the total separation is 
a minimum. Chains of corrective operations connect- 
ing matched pairs can then be applied. Prior work 
has calculated the separation of two syndrome changes 
si = s 2 = fejiiA) using d(si,s 2 ) = \i\ - 
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\ti — t<i I Tfj, [13|, [14( , however it was shown 



that this is far from optimal and leads to poor per- 
formance, particularly at low error rates. In this work, we 
instead approximate the probability P(s\,S2) of a given 
pair of syndrome changes being connected by an error 
chain, and set d{&\,&%) = — ln(P(si, Sa))- This choice of 
rf(si, S2) is natural, accounting for the substantial perfor- 
mance increase we observe. 

To calculate P(si , S2), we must study the effect gate er- 
rors. Fig. [3] shows all possible pairs of syndrome changes 
resulting from all possible errors on all meaningfully dis- 
tinct gates. The CNOTs shown measure an X-stabilizcr. 
The effect of errors on the CNOTs used to measure a 
Z-stabilizer can be obtained by interchanging X and Z . 

Using fig. [3j fig. [4] was constructed, grouping gate er- 
rors leading to specific pairs of syndrome changes. The 
probability of an odd number of errors occurring in each 
group gives the probability of the associated link. Us- 
ing the standard error model, the probability of the link 
shown in fig. [4jt is 
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Defining similar probabilities ps, PCi Pd, Pe, Pf for 
fig. Hfc-f, fig. [5] shows the 0(p) links from one syndrome 
to its neighbors. Some straightforward modifications of 
the links and expressions are required at the temporal 
and spatial boundaries. 

The probability P(si,S2) that two syndrome changes 
are connected is the sum of the probabilities of all con- 
necting paths. The probability of a given path is the 
product of the link probabilities along the path. Several 
approximations of P(si,s 2 ) are worthy of study. The 
simplest approximation is to take a single path of max- 
imum probability P m ax(si,S2) and define dmax(si, S2) = 
— ln(P max (si, S2)). We shall see that this approximation 
is sufficient to substantially increase pth, and that more 
accurate approximations do not lead to further increase. 
The performance of surface code QEC using d ma x(si, s 2 ) 
and the standard error model is shown in figs. [SH7] 
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FIG. 3: Syndrome changes resulting from a-d) specific two- 
qubit errors on specific CNOTs, e) a syndrome qubit mea- 
surement error, f) a data qubit memory error. 
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FIG. 4: Numbered error processes from fig. [3] contributing to 
specific links. Superscripts 1, 2 and 3 indicate errors occur- 
ring with probability 8p2/15, 2p//3 and pm respectively. All 
others occur with probability 4p2/15. 



FIG. 7: Average rounds of error correction before logical 
Z failure as a function of the gate error rate p when using 
dmax(si, S2) and the standard error model. 




FIG. 5: All possible links from a syndrome change to its neigh- 
bors. Letters correspond to figs. UJi-f. 



Figs. |HH3 give strong evidence of p t h = 1.1%. We 
have verified this by simulating high distance codes at 
p = 1.1% and observing neither increase nor decrease of 
the failure time. This is enormously encouraging, and 
motivates one to better approximate P(si,S2) in an ef- 
fort to further increase pth- Additional accuracy can be 
achieved by taking all shortest length paths (measured in 
links) between s± and S2 and calculating the sum of prod- 
ucts of link probabilities along each path. We shall define 
the resulting distance measure as do(si,S2). We can de- 
fine similar distance measures d n (si,82) taking into ac- 
count all minimum length / paths and paths of length no 
greater than l + n. The performance of surface code QEC 
around p t h using do( s i ; ^2) is shown in fig. [8] It can be 
seen that p t h remains 1.1%. 

The fact that do results in the same pth as using a sin- 
gle maximum probability path distance measure c? max can 
be explained by noting that do(si,Sa) oni y differs from 
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FIG. 8: Average rounds of error correction before logical 
X failure as a function of the gate error rate p when using 
do(si, S2) and the standard error model. p t h remains 1.1%. 



^max(si, S2) if si, S2 are separated by at least two links. 
Single link paths are unique minimum length paths im- 
plying do = d max - The vast majority of error chains, even 
for p = pth, are single links. We find that the modifica- 
tion of the distance associated with the multiple paths of 
syndrome changes separated by multiple links is in the 
region of 10-20%. Given such multiple link paths are 
not leading order contributors to pth in the first place, 
this relatively small weight change does not result in an 
observable improvement of pth- 

Higher order approximations d n (si, S2) will also result 
in the same pth as the distance is hardly altered by in- 
creasing n. To take a typical example, for p = 0.01, 
si = (0,0,0), s 2 = (1,-1,0) we obtain d = 6.91 (6 
paths), d\ — 6.86 (30 paths) and d^ = 6.85 (390 paths). 
The exponential increase of the number of paths is well 
balanced by the exponential decrease of the probability 
of these paths. 

The balanced error model is a better model of all quan- 
tum gates failing with equal probability than the stan- 
dard error model, and appropriate modification of the 
polynomials using P2 = p, pi = 4p/5 and pu = 8p/15 
leads to p t h = 1.2%. The ion trap error model, with p 2 = 
p, pi = p/1000 and pm = p/100 leads to p t h = 1.4%. 
Arbitrary stochastic error models are straightforward to 
analyze using our formalism. 

To conclude, by performing a detailed study of the 
probability of different pairs of syndrome changes and 
feeding the simplest approximation of this information 
into the minimum weight perfect matching algorithm, 
we have been able to raise the geometrically constrained 
threshold error rate to 1.1-1.4%, depending on the exact 
error model, while maintaining computational efficiency. 
This is the first time a geometrically constrained thresh- 
old error rate has been observed over the 1% level. There 
is the potential for still further improvement by taking 
into account correlations between X and Z errors, which 



shall be pursued in further work. 
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